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SNYTHESIS OF REGIONAL CRUST AND UPPER-MANTLE STRUCTURE 
FROM SEISMIC AND GRAVITY DATA 

Abstract 

The principal objective of this study is to combine available 
seismic and ground-based gravity data to infer the three-dimensional 
crust and upper mantle structure in selected regions. This snythesis 
and interpretation proceeds from large-scale average models suitable 
for early comparison with high-altitude satellite potential field data 
to more detaile.’ delineation of structural boundaries and other vari- 
ations that may be significant in natural resource assessment. While 
the study focuses primarily on seismic and ground based gravity data, 
other relevant information (e.g. magnetic field, heat flow, Landsat 
imagery, geodetic leveling, natural resources maps) ure used to constrain 
the structure inferred and to assist in defining structural domains and 
boundaries. The relevant seismic data base consists of regional 
refraction lines, limited reflection coverage, surface wave dispersion, 
teleseismic P and S-wave delay times, anelastic absorption (Q) , and 
regional seismicity patterns. The gravity data base consists of 
available point gravity determinations for the areas considered. The 
primary area considered is the Eastern United States from the 
Mississippi River to the Atlantic Continental Margin. 

The interpretation makes use of modern inversion methods, digital 


analysis techniques, and empirical evidence on density- seismic velocity 
relationships for crustal rocks. Results of the first year's study 
under NASA Grant No. NSG 527b are summarized in this report. 



SUMMARY 


The first year of this investigation (under Grant NSG 5276) has 
been devoted heavily to reviewing the available data sources, gathering 
previous regional seismic results and commencing a snythesis of gravity 
(free air and Bouguer) anomaly data and seismic refraction and reflection 
information. (The second year of this research project originally 
proposed for a three-year interval is being continued under NASA Grant 
NCC 5-19.) 

To look for possible large-scale intraplate crustal features in 
eastern North America that may be associated with past or current tec- 
tonic structural elements, we analyzed 1° x 1° free-air anomaly data. 

This analysis included conventional and new techniques developed in 
this study for enhancement of structural boundaries. Several large, 
regional blocks (> 500 km in lateral dimension) were identified that 
likely reflect deep crustal or upper mantle density differences. 

However, the correlation of regional seismicity with these large block 
structures is not strong enough to be compelling, although the major 
gravity anomaly patterns appear to follow major tectonic trends in the 
Eastern United States. Several smaller, distinct structural blocks are 
likely contained within each large block and these may be important to 
the tectonic interpretation. The details of this investigation are given 
in the Appendix (Abraham Biadgelgne's M.S. paper). 

We have also obtained and reviewed especially carefully the data 
from a comprehensive tectonic study of the Southern Appalachian Region 
recently conducted by the Tennessee Valley Authority (TVA) that includes 
gravity, seismic, remote sensing, aeromagnetic , and regional geologic 


observations. These and other dutu that we have studied in detail 
suggest the presence of distinct crustal blocks that have lateral 
dimensions of the order of 50 km, and we are attempting to determine 
whether this characterization is typical of most ureas in the EUS. 
Contained in the TVA study area is a major magnetic source body in 
S. Central Kentucky with associated strong gravity and seismic anomalies 
that extend upward to the base of the Paleozoic section; this feature 
and others that we can identify such as the Scranton gravity high are 
receiving special attention. In particular these major crustal anomalies 
are being investigated and modeled in more detail by K. Hawman in a 
M.S. thesis that is in progress. Details of his investigation will be 
reported under the continuation of this study (NASA Grant NCC 5-19). 

During the past year we have also communicated frequently with 
NASA (Goddard) scientists to coordinate our work with the MAGSAT effort 
and we plan to continue this cooperative exchange. A vist to Penn State 
by Drs. Langel , Taylor, Allenby and Thomas was especially helpful in 
this regard. 

Based on this initial year's study, future efforts should: 

(1) concentrate on obtaining more detailed Information on the crustal 
structure in those areas exhibiting .zrjor (MAGSAT) regional magnetic 
anomalies, (2) refine models of the crustal structure in areas where 
seismic results (interpretations) are now lacking, and (3) combine our 
models with MAGSAT results. Crustal structure estimates from regional 
seismic networks that are now being operated in the EUS should be 


incorporated as well. 
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ABSTRACT 


An interpretatlonal technique baaed on the third 
vertical derivative of a gravity anomaly has been developed. 
The versatility and possible limitations of the method were 
evaluated by examining the response for two- and three- 
dimensional rectangular prisms with a wide range of sizes 
and depths. The effect of the Interference of adjacent 
sources was also studied. It was demonstrated that for 
Isolated prismatic models with thicknesses ranging from 0.2 
to 20.0 sampling Intervals, the third vertical derivative 
of the g. ivity anomaly outlines the top surface of the 
sources with error bounds of less than 8Z. For cases In 
which the sources were in contact, the boundaries were 
delineated by well defined low or high trends. Moreover, 
the horizontal distance between the maxima and minima of the 
vertical derivatives gives an estimate of the depth to the 
source with error bounds of less than 9X If the thickness Is 
less than two grid intervals. Depth estimation for sources 
In contact was complicated by the interference of the 
extrema of the derivatives. An approach for Improving the 
results In such cases was proposed. 

An attempt was then made to enhance some linear 
features In the most recent 1° by 1° gravity data for the 


central and eastern United States using the third vertical 



Hi 


derivative 
lineaments - 
assisted in 
problem of 


technique. Correlation of the Inferred 
1th some seismic trends and magnetic features 
outlxng soma likely crjstal blocks. The 
the cuperpositlon of the effects of regional 


shallow and deep structures was discussed. 
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CHAPTER I 


INTRODUCTION 

Quantitative Interpretation of data collected in a 
gravity or magnetic field survey involves investigating 
various types of models that will best fit the observed 
anomaly and the known geology of the region of interest* 
For this reason, there is no paucity of papers discussing 
techniques for evaluating the field due to different model 
sources (e.g., Henderson and Zeitz, 1948; Vacquler, 1958; 
Battacharyya , 1964, ]96f). From among such models, the 
rectangular block-shaped body has found the most widespread 
applications in the Interpretation of potential field data. 
The popularity of the model arises from its versatility in 
representing a wide variety of geological structures and 
from its utility as a building block for complex mass 
distributions. Looking into the future, one expects the 
model to gain more importance as satellite determination of 
the earth's gravity and magnetic field Increases the 
desirability and the feasibility of studying major crustal 
structures. 

In the past, various qualitative and quantitative 
interpretat lonal techniques have been developed to assist in 
determining the physical parameters of field sources. Of 
all these, a method that has been used extensively to 
interpret aeromagnetlc data from various geological 
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provinces is the one first developed by Vacquier et al • 
(1958) and later extended by Battacharyya (1964, 1966) and 
Jain, et al . (1974). Vacquier, et al (1958) conducted 
detailed model studies and showed that the zero contour of 
the second vertical derivative of the magnetic field closely 
outlines the top surface of a vertically sided, isolated 
prismatic model with infinite thickness and vertical 
magnetization. Moreover, Jain, et al . ( 1974 ) shoved that 
the horizontal distance between the plotted maxima and 
minima of the derivatives gives a reasonably accurate 
estimate of the depth to the source. 

Standard lnterpretat ional techniques have been applied 
for the gravity case. However, in the past, no attempt has 
been made to apply the Vacquier, et al . ( 1958) approach. 
Hence, in this paper it will be shown that the general 
procedure of determining the horizontal dimensions and the 
depth of a prismatic source is also applicable to gravity, 
since Poisson's relation implies that the second vertical 
derivative of the magnetic field of a source is equivalent 
to the third vertical derivative of its gravity field, 
provided certain basic assumptions discussed later are met. 
Two major differences between the gravity and the magnetic 
problems are expected to affect the results of the gravity 
case. In magnetics, the prismatic models are assumed to 
have infinite thicknesses. This allows one to deal only 
with a sheet of monopoles at the top surface. However, in 
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gravity, there Is a continuous distribution of density with 
depth which requires the model be limited to a finite 
thickness In order to avoid a singularity* Thus, the effect 
of finite thickness on the position of the zero contour and 
on the estimation of depth must be determined* An Important 
advantage of using gravity data over magnetic data Is that 
there Is no need to correct field data for inclined field 
direction, l.e*, the gravity field is already "reduced to 
the pole." The Implication of these basic factors will be 
explored In detail In later sections for both two- and 
three-dimensional model'* For the sake of simplicity, most 
of the analysis will be conducted In the wavenumber domain* 
Finally, the interpret at lonal technique will be applied 
to the most recent 1° by 1° free-alr gravity data of the 
eastern and central United States prepared by Rapp(I978). 
The main Intent of this study will be to enhance some major 
linear features suggested by the field data* Some basic 
limitations of the data set. that make the final outcome of 
this aspect of the investigation preliminary in nature will 


also be addressed 



CHAPTER II 


THEORY 

The theoretical framework for the lnterpretat lonal 
technique developed in this paper will be presented herein. 
The general approach will be to begin with a review of the 
lnterpretat lonal scheme used In magnetics, and then to 
invoke Poisson's theorem to show the nature of the operation 
for gravity data* This will be followed by a discussion on 
the gravity effects of a rectangular prismatic model with 
finite thickness. An analysis of both the continuous and 
discrete forms of the vertical derivative operators will 
also be presented in the wavenumber domain for later use. 

Magnet ic Theory and Its Applications 

Vacquier, et al . (1958) have shown that for bottomless 
rectangular prisms that have great oreal extent compared to 
their depth of burial, the zero contour for the second 
vertical derivative of the magnetic field closely follows 
the outlines of the top surface of an isolated body. The 
theoretical framework for this concept is based on the 
assumption that the anomalous field results solely from 
induction in the earth's magnetic field and that the field 
is nearly vertical. However, experimental work on rock 
magnetism by several investigators (Green, 1960; Ross and 
Lavln, 1966; and Bat tacharyya , 1964 ) has demonstrated that 
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remanent magnetization Is commonly present and may even 
dominate the Induced magnetization* Despite this 

restriction, a wide application of this rather simple 
technique over diverse geological provinces has lead to 
significant results (Jain et al . , 1974 ; Bhat tacharyya , 1966; 
Zeltz et al . , 1967)* In cases where a relatively strong 

remanent magnetization exists at an angle to the earth's 
field, or where the dip of the magnetization vector Is too 
low, the contours would be skewed to one end of the model 
and the Inflection points would not mark the edges of the 
model (Zeltz and Andreasen, 1967)* For such cases, the 
asymmetry in the contours may be improved by reduction to 
the pole if the magnetization direction is known. 

The other main information obtained from the vertical" 

derivative operator, in magnetics, is the estimation of the 

depth to the top of the prismatic source. The horizontal 

distance between the maxima and minima of the second 

vertical derivative, R, has been related to the depth of the 

source, D , through theoretical studies and detailed 
s 

empirical observations. Jain et al . (1974) have shown that 

R varies in almost direct proportion to D g , or 

D - CR ( 1 ) 

s 

where C is a constant that ranges from 1.2 for very narrow 
or very wide bodies to about 1.4 for bodies whose width is 
double the depth. In general, taking the vclue 1.3 for C 
limits the percentage errors in depth estimation to within 
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1 0 Z > Once again. If the reduction to the pole Is 
Incomplete, due to a lack of information on the direction of 
magnetization or due to inadequate data processing, the 
value of C may be as high as 1.5 leading to erroneous depth 
determinations (Jain, 1974). It has also been observed that 
for cases in which the distance between different sources is 
roughly equal to or less than the depth of burial, the 
technique provides less reliable depth estimates. 

Poisson's Theorem and Its Appl lcatlons 

The Newtonian or gravitational potential and the 
magnetic dipolar potential are related by a simple equation 
provided the restrictive assumptions discussed below are 
met. Suppose the mass of a body that fills a volume, V, has 
a continuously distributed dipole moment per unit volume, 
I?, where I is its magnitude and "t is a unit vector in the 
direction of magnetization. Then, the magnetic potential, 
A( r) , at an arbitrary point, S, outside the body (Figure 1A) 
is given by 



Similarly, the gravitational potential due to a continuous 
mass of density, D, at the same exterior point, S, is given 
by 

U(f) - - QD J , (3) 




A. Spatial Domain 



B. Wavenumber Domain 


Figure 1. Coordinate System. 
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•8 

where Q is the Gravitational Constant (6.67 x 10 cgs 
units). Assuming that the density, D, is uniform and the 
magnetization, I, is uniform both in magnitude and direction 
over the same volume, V, tie geometric factor of Equations 2 
and 3 may be equated giving 

A< *> • q"d If u<?) ■ (4) 

The above equation, which is attributed to Poisson (1826), 
implies that all the properties of the magnetic potential of 
an anomalous body are derivable from its gravitational 
potential or vice versa, provided the basic assumptions are 
met. For cases in which the vertical component of the 
magnetic field are considered, Equation 4 ma"y be 
differentiated to obtain 

H - 1 9 

2 Q D 3t 8 z ’ (5) 

where and g^ are the vertical components of the magnetic 
and gravity fields, respectively. 

Equations 4 and 5 have been widely applied in the 
analysis of local and regional anomalies. Garland (195ib) 
used the relation in transforming the Crow Lake gravity 
anomaly in the Canadian Shield into a magnetic anomaly, and 
demonstrated that the source was non-unif ormily magnetized. 
Lundback (1956), Ross and Lavin (1966), and Robinson (1971) 
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have applied various versions of Poisson's relation to 
determine the direction of magnetization for two- and three- 
dimensional bodies by the transformation of gravity or 
magnetic fields* Recently, a numerical approach that 
employs least-square inversion of gravity and magnetic data 
in the wave number domain has been successfully used to 
determine the ratio I/D, and the direction of magnetization 
(Cordell and Taylor, 1971)* Kanasewich and Agarwal (1970) 
have also applied the theorem in the wavenumber domain to 
compare gravity data measured over a large area with the 
corresponding magnetic field data reduced to the pole* 

In the present study, Poisson's relation is used to 
relate higher order vertical derivatives of gravity and 
magnetic fields* Suppose that the field source is 
vertically magnetized, or that it is reduced to the pole. 
Then, Equation 5 becomes 


H 

z 


• I _ 3 _ 

Q D 3z %z 


( 6 ) 


Differentiating twice with respect to z gives 
a 2 _ I 9’ 

3z 2 2 Q D 3 Z 3 8 z 


Equation 7 shows the second vertical derivative of a 
magnetic field corresponds to the third vertical derivative 
of the gravity field* Thus, one can extend the magnetic 
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analysis procedure of Vacquier, et al. (1958) and Jain, et 
al . (1974) fo gravity anomalies by employing the third 
vertical derivative of the field* 


Gravltat lonal Anomaly o f Rectangul ar Priam 

A right rectangular prlam with sides parallel to the 
coordinate axis, at depth H below the surface of observation 
is shown in Figure 2* The derivation of a closed expression 
for the vertical component of the gravitation attraction at 
an arbitrary po in t , S (x , y , z ) , outside of or on the boundary 
of the prism was given by Nagy (1966)* Through a simple 
transformation of his coordinate system, the general 
equation for the field value at an arbitrary point becomes 


g 

z 


QD 


(a - x)ln(g - y + R) + (6 - y)| n (a - x + R) 


- (Y - z) arcs in 


2 2 2 

(6 - y) + (y - *) + (b - y) R 

(6 - y + R) (6- y) 2 + (y - z) 2 


(S) 


where the variables are evaluated from X - A to X + A for 

c c 

at, from Y - B to Y + B for A and from H to H + T for 1 
c cr 

(Figure 2). 

The complexity of the space-domain expression limits 
its usefulness in in te rpr e t a t iona 1 schemes for all practical 
purposes* All the parameters of the prism are contained in 
coupled terms and extraction of any of them appears rather 
difficult* Moreover, the analytic evaluation of the third 
vertical derivative in the spatial domain would be very 
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cumbersome at best. The first vertical derivative alone 
Involves over tea terms, end it is not difficult to Imagine 
how Intractable the third vertical derivative would be. 

Alternatively, the gravity effect may be represented by 
its Fourier transform given by 


C(u,v) 


j"g(x,y) exp(-i(ux + vy) dx dy (9) 

•y 


where 


g(x,y) ■ 'A IT' 1 ! j G(u,v) exp(iux + ivy) dx dy , 


(] 0 ) 


and u and v are wavenumbers in radians per unit length in 
the X and Y directions, respectively. 

It is worthwhile to note that to have a definite 
Fourier transform, in strict mathematical sense, a function 
needs to be absolutely lntegrable over the X-Y plane. This 
requires the function to go to zero over large values of X 
and Y. However, the gravity field of anomalous sources does 
not vanish over the surface of observation, and, lu fact, 
theoretically, it extends to infinity. To overcome this 
mathematical difficulty, one must assume that ihe field will 
be defined over 3 large, but finite area of the coordinate 
plane and that it becomes zero outside this domain. A 
method for reducing truncation effects that arise from such 
finiteness of the surface of observation will be addressed 
in the data processing section of the next chapter. 
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Carroll (1973) showed that the spectrum for a 
rectangular prismatic model with finite thickness la given 
by 

G( u v ) - 2TfC<M exp(-Hp) (1 - exp(-Tp) ) sinc(Au) 

T p 

sinc(Bv) exr(-ix c u) exp(-iy c v) . (11) 

2 2 2 

where M Is the total mass and p ■ u + v (Figure IB). The 
wavenumber expression of the gravity field is simpler to 
analyze since the parameters are now uncoupled Into separate 
terms. The term exp(-Hp) represents the depth term, (1 - 
exp(-Tp))/Tp represents the thickness term, and the sine 
functions, slnc(Au) and slnc(Bv) represent the width terms. 
It is also noted that Equation 11 may be modified further to 
accomodate two-dimensional problems, by assuming the gravity 
field to be constant along one of the axes. If, for 

Instance, the field is constant parallel to the Y axis, the 
v dependence in Equation 11 will drop out, leaving the 
spectrum for a rectangular pc r al 1 ep 1 ped extended along the Y 
axis to infinity (Carroll, 1973), or 

2ttQM 2 

G(U,V> “ ~ x | U | ex P<-H|u|)(l - exp(-T |u| ) ) 

jsinc(Au) | 6 (v) , (12) 

where M 2 la the mass per unit length of the source and S(v> 
is the delta function. 
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Discrete Form of the Nth Vert leal Derlvat lve 

The vertical derivative of order n has a wavenumber 

response given by (Agarwal, 1972) 

F n (u, v) - (u 2 +v 2 ) n/2 • (13) 

In particular, for our purpose, the ideal third 

vertical derivative operator has the form 

F 3 (u, v) - (u 2 +v 2 ) 3/2 . (14) 

Hence, the third vertical derivative of a gravity field 

3 

can be easily obtained by multiplying G(u, v) by F (u, v) as 
shown below 

G 3 ( u , v) * G(u, v) F-(u, v) , (15) 

and then inverse Fourier transforming the result. The 
simplicity of the approach in the wavenumber domain stands 
in clear contrast to the rather intractable alternative of 
the analytical evaluation of the derivative in the spatial 
domain . 

The gravity and magnetic effects that observed field 

data represent are, in reality, continuous functions, while, 

in practice, one often only knows these functions at points 

of intersection of a grid. This requires that we discretize 

the smooth and continuous vertical derivative operator, 
3 

F (u, v). If we have N data points in the X-direction and M 
data points in the Y-direction, the sampling intervals in 
the wavenumber domain, in cycles per data interval will be 


] /n and Ak 2 = 1/ M 


(16) 


15 


Equation 14 then becomes 

oVpkJ - (( -2Sl ) 2 + ( _2Hk 2 ) 2 ) 3/2 U7) 

1 1 N M * 

where k^ and k^ are the Integers 0, 1, 2 For the 

discrete case Equation 15 may be rewritten as 

G 3 (k 1 ,k 2 ) - G(k lt k 2 ) D 3 (k lt k 2 ) . (18) 

It Is Important to note that the derivatives (Equation 
17) are Independent of the Internal mass distribution and 
can be computed unambiguously for any observed surface data. 
Moreover, the form of the vertical derivative operator shows 
the tendency to amplify high wavenumber anomalies and 
observational errors relative to broad features In the 
surface field. It Is, thus. Imperative to make a detailed 
Inspection of the spectrum of the measured anomaly, and then 
decide which wavenumber range should be operated on and 
which parts of the Fourier spectrum should be suppressed, 
thereby gaining full control over the output of the 


derivative operator 
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Determination of W idth and Depth Parameters 

In the previous sections we have shown that Poisson's 
theorem justifies an approach for gravity field analysis 
equivalent to the magnetic lnt e rp re t a t lonal technique. This 
implies that the zero contour of the third vertical 
derivative of a gravity field will closely outline the top 
surface of a rectangular prism with finite thickness. In 
cases where a continuous distribution of sources is 
encountered, the zero contours would be expected to 
delineate the boundaries between the sources. In fact, one 
apparent advantage In gravity Interpretation is that the 
anomaly Is "naturally" reduced to the pole. Hence, there is 
no problem of the zero-contour being skewed to one edge of 
the prismatic models, as would be the case In magnetics when 
the source Is not vertically magnetized. On the other hand, 
care must be taken in specifying the thickness of the 
models. In magnetics, the effect of the north and south 
poles of the magnetization vectors of a vertically-sided 
prismatic model cancel each other except for the ones at its 
top and bottom surfaces. This implies that the vertical 
component of the magnetic intensity at any point of 
observation is proportional to the difference of the solid 
angles subtended by the top and bottom surfaces, 
respectively. Hence, in the case of bottomJess rectangular 
blocks, one has to consider only the effect of a sheet of 
monopoles at the top. However, in gravity, there is a 
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continuous distribution of density with depth, and extending 
the thickness to Infinity will lead to a singularity. This 
requires that, unlike the approaches followed in magnetics 
(Vacquier, et al., 1958; Brattacheryya , 1964, 1966;Jain, 
1974), the gravity investigation must be restricted to 
prismatic models with finite thicknesses. Hence, one must 
study the relation of thickness to the position of the zero 
contour and determine the range of thicknesses over which 
the method works. 

Another major application of the vertical derivative 
scheme is based on the relationship of the depth to the top 
of the source and the horizontal distance between the maxima 
and minima of the second vertical derivative of the magnetic 
field as shown in Equation 1 (Jain, et al., 1974"). It 
should also be noted that Equation 7 gives a linear 
relationship between the second vertical derivative of a 
magnetic field and the third vertical derivative of a 
corresponding gravity field. Hence, the method for 
estimating the depth to the source is readily transferable 
to the interpretation of an anomaly due to an isolated 
gravity source, provided the implications of the 
restrictions put on the thickness as discussed above are 
respected. 


CHAPTER III 


PROCESSING AND INTERPRETATION ON MODEL DATA 

To examine the utility of the lnterpretat lonal 
technique formulated In Chapter II, the gravity anomalies of 
several numerical models and their third vertical 
derivatives were studied using various data processing 
schemes. For the sake of simplicity and ease of study, most 
of the detailed analysis was focused on two-dimensional 
rectangular parallelepipeds with finite thicknesses. This 
was supplemented by a limited consideration of three- 
dimensional models, to show the effect of having two or more 
models that are in contact or well Isolated. The overall 
reliability and possible limitations of the method were also 
evaluated. Several of the more Interesting two- and three- 
dimensional model studies are described herein. 

Data Processing 

The gravity anomaly for a right rectangular prism of 
finite thickness was generated using computer programs based 
on the Talawanl algorithms (Talawanl, 1965; Talawanl and 
Heirtzler, 1964). The anomaly was then f a 3 c- Four 1 e r 
transformed (Cooley and Tukey, 1965) for processing with the 
discrete form of the third vertical derivative operator 
(Equation 17). 
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A problem was encountered early In the study brought 
about by the basic assumption that the field becomes zero 
outside the area of observation* Such data-set truncation 
In the spatial domain results In high wavenumber 
oscillations In the wavenumber domain. To avoid this 
undesirable phenomenon, a Hanning window was found to be the 
most suitable pre-tr ansf orm filter. The form of the window 
used for a two-dimensional data set Is 

W(x,y) ■ -r~ (1 + cos ttx) (1 + cos Try) 

4 M” M (19) 

for x <N and y < M , 

where N and M are the array sizes of the gravity field data 
In the X and Y direction, respectively. For profile data, 
one of the variables (X or Y) is set to zero. 

Another problem encountered early in the project was 
the sensitivity of the vertical derivative operator to the 
selection of grid spacing. Model profiles were generated 
with even spatial sampling Intervals of different lengths. 
In general, the results demonstrate that the finer the grid 
spacing, the more accurate the computed values of the 
horizontal dimension and the depth of the prismatic models. 
On the other hand, It Is important to note that the 
superimposed grid acts as a filter that retains anomalies 
pertaining to sources of certain wavelengths. However, in 
some geophysical studies local or shallow anomalies are of 
minor or no Interest and should be eliminated. Hence, in 
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practice, there Is a need to use a grid spacing which Is 
commensurate with the structural features being studied 
provided alaislng of the data-set Is avoided. For our 
particular purpose, a sampling Interval of about 0.2 times 
the width of the numerical modelu gave adequate results In 
outlining the top surface of the sources , while a value of 
0.1 was required for estimating their depth of burial. 
These values result In oversampled data. In cases where 
coarser grid spaclngs were employed. Interpolation to finer 
spacing using the bi-cublc spline interpolation In the 
spatial domain and padding with zeroes In the wavenumber 
domain were found to be particularly useful In improving the 
results. 

Delineation of Block Boundaries 

The effectiveness of the zero contour of the third 
vertical derivative in outlining the edges of an isolated 
source or In delineating the lithologic boundaries of 
continuously distributed sources are influenced by the 
physical parameters of the bodies. Detailed and systematic 
empirical Investigations were carried out to determine 
possible restrictions introduced by these factors. Some of 
the most signiflca.it results will be discussed below. 

One of the prerequisites for carrying out a field 
survey or a model study is selecting an appropriate sampling 
scheme. In field applications, the grid spacing would be 
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set depending on the data coverage desired and/or the 
logistics of the date acquisition process. In model 
studies, they are selected to give reasonably accurate 
results without being too fine for practical applications. 
The proper grid interval selected on this basis has been 
discussed in the data processing section. In the subsequent 
sections, the procedure followed is to use an adequate 

sampling interval (0.2 or 0.1 times the width) for a 
particular model and then to vary the physical 
pa rame te r s , keeping the grid spacing fixed. In this manner, 
one determines the range of parameters that will be covered 
for a given condition without introducing serious errors. 
The standard model used for this purpose has a width equal 
to 12.0 sampling intervals. The parameters were varied over 
ranges consistent with the nature of the gravity field data 
expected in the next chapter. It should also be noted that, 
from now on, the grid spacing is the unit of measurement 
when it is not specified otherwise. 

Two-Dimensional Models 

Figure 3 illustrates a typical gravity anomaly and its 

third vertical derivative for a two-dimensional rectangular 

paralle lplped . The location of the maxima and minima of the 

0 0 

derivative are indicated by A and B , respectively. The 

position of the reversal of sign in the values of the 

0 

derivative is indicated by C . This corresponds to the zero 
contour in the three-dimensional case. 


MGAL o , MGAL/DI 
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. Third vertical derivative of the field 



B. Gravity anomaly 


A. Model geometry. 


Figure 3. The gravity anomaly and third vertical derivative 
for a rectangular prism. 
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Table 1 gives a sample of interpretational results 
showing the effects of thickness on width estimations. 

A two-dimensional rectangular paralle lpiped with a width 
equal to 12.0 was employed. The depth was kept fixed at 
2.0, while the thickness varied from 0.2 to 20.0. The data 
indicates that accurate determinations of width with error 
bounds of less than 2.2Z are to be expected over all the 
range of thicknesses considered. A slight increase in the 
differences between actual and predicted width is noticed as 
the thickness increases by two orders of magnitude. It 
should be noted that, since there is a continuous 
distribution of density with depth, increasing the thickness 
amounts to adding some more mass to its deep end. Hence, 
in the limit as the thickness goes to infinity one would be 
faced with a singularity. This would, necessarily, lead to 
a slowly increasing inaccuracy in the results. 

Table 2 shows the effect of using a fixed grid interval 
and profile length that were found adequate for a model 
width equal to 12.0 for models of other widths. 

The depth for all models was 2.0, and the range of width 
considered varied from 2.0 to 24.0. The thickness was kept 
constant at 4.0. Truncation effects introduced as the 
models extend closer to the border probably contributed to 
the slightly higher errors for wider models, though the 
Hanning window has greatly reduced such effects. The 
narrower models, on the other hand, have been greatly 
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TABLE 1 

EFFECT OF THICKNESS ON WIDTH ESTIMATION 


Es tinated 

Thickness Width Percentage 

( sampling unit ) ( sampling unit ) Error 


0.2 

11.88 

-0.7 

0.4 

11.86 

-1.2 

0.6 

11.84 

-1.3 

0.8 

11.84 

-1 . 3 

1.0 

11.82 

-1.5 

1.2 

11.82 

-1.5 

1.6 

11.80 

-1.6 

2.0 

11.80 

-1 .* 

4.0 

11.76 

-2.0 

8.0 

11.74 

-2.2 

12.0 

11.74 

-2.2 

16.0 

11.74 

-2.2 

20.0 

11.74 

-2.2 


affected. The smallest model presented shows an error of 
30%. This primarily reflects the effect of inadequate 
sampling of the data for the particular case. The same 
analysis was repeated for a range of thicknesses extending 
from 0.2 to 20.0 and a similar trend in percentage errors 
was observed. The only difference was that the thickness 
effect makes the overall errors slightly higher as was the 
case in Table 1. 

Several models were evaluated to examine the influence of 
the depth parameter in outlining the top surface (Table 3). 

Within the range of practical interest for local and 
regional studies, the results are reasonably accurate. As 
we probe for deeper sources, the inaccuracy in the predicted 
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TABLE 2 


WIDTH ESTIMATIONS 

FOR FIXED SAMPLING 

INTERVAL 


Es tlmated 


Width 

Width 

Percentage 

(sampling unit) 

(sampling unit) 

Error 

2.0 

2.60 

• +30.0* 

4.0 

4.24 

+ 6.0 

6.0 

5.86 

- 2.3 

12.0 

11.76 

- 2.0 

18.0 

17.60 

- 2.2 

24.0 

23.42 

- 2.4 


* curve very oscillatory 

values increases slightly. In fact, the errors become even 
larger when the effect of width .and thickness shown in 
Tables 1 and 2 are included with the depth effect. This 
would be especially true for very narrow sources that are 
extremely thick. Hence, once again, it is noted that the 
extremely narrow sources suffer more than the extended ones. 
This is, perhaps, reflective of the nature of the depth term 
in Equation 12, together with the problem of Inadequate 
sampling discussed before. The data spectrum in Equation 12 
falls off as exp(-H|u| ) with increasing value of H. Thus, 
the depth factor acts as a low-pass filter. 

It should be noted that the results given in the three 
tables cover a wide range of depth, width and thickness that 
are quite realistic. If, for instance, a grid interval in a 
field survey is assumed to be 10 km, the range of 
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TABLE ^ 

EFFECT OF DEFTH ON WIDTH ESTIMATION 
Estimated 

Depth Width Percetage 

( samp ling unit ) ( sampling unit ) Error 


1.0 

• 11.94* 

-0.52 

2.0 

11.80 

-1.72 

3.0 

11.60 

-3.32 

4.0 

11.44 

-3.32 

5.0 

11.26 

-6.22 

6.0 

11.10* 

-7.52 

7.0 

11.04** 

-8.02 

curve very 

oscillatory 



thicknesses (Table 1) ranges from 2 km to 200 km. The range 
of widths extends from 20 km to 240 km, while the depth 
varies from 10 km to 70 km. Similar range estimations are 
plausible for other grid intervals depending on the field 
survey desired. It is also noted that for most of the 
ncise-f ree data that has been considered, the range of 
errors is less than 8 percent. These error values Indicate 
that the results are reasonably accurate for practical 
purposes . 

Three-Dimensional Mode Is 

To illustrate the versatility of the int e rpre tat Iona 1 
technique in outlining Isolated gravity sources, and in 
delineating lithologic boundaries, several empirical models 
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were run using a three-dimensional lactangular prism with 
finite thickness as the basic model. Some of the more 
interesting results are discussed below. 

The effect of truncation was studied by processing a 
segment of the anomaly generated for a rectangular prism of 
finite thickness (Figure 4). Several pre-transform filters 
were tested in an attempt to smooth the data spectrum and 
prevent spectral leakage. However, the output of the 
vertical derivative operator persistently showed some 
instability at the truncated end. To circumvent the 
problem, the derivative operator was designed in the spatial 
domain. The set of weights were derived by inver se-Fourier 
transforming the wavenumber expression of the third vertical 
derivative operator and then truncating it into a set of 7 
by 7 coefficients by using a symmetric Hanning window 
(Equation 1?). The wavenumber response of the filter is 
shown in Figure 5. It is noted that over a wide range of 
wavenumbers the behavior of the amplitude response is 
circularly symmetric as it would be with the ideal case. A 
reversal of the curvature of high wavenumber contours occurs 
for wavenumbers greater than about C.4 cycles per data 
interval (cycles/di). This effect, which is also observed 
in several of the wavenumber responses of operators reported 
in the geophysical literature (Henderson, 1960; Elkins, 
1951; Fuller, 1967) is primarily a result of the loss of 
circular symmetry caused by specifying the discrete weights 
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Figure 4. A segment of the gravity anomaly generated for a 
three-dimensional rectangular prism (contour intervals in mgals). 
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Figure 5. The amplitude response of the third vertical 

derivative operator designed in the spatial domain, 
(contours interval = 3 db ) 
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on a square grid (Fuller, 1967). However, most of the 
energy content of the data is concentrated below about 0.1 
cycles/di, and should not be affected by the asymmetry. 
Considerable improvement in the accuracy of the output was 
achieved by convolving the set of coefficients with the 
gravity anomaly (Figure 6). The resolution of the upper 
corner of the body and the delineation of the truncated end 
is rather striking. It should also be noted that the loss 
of data around the edges during the convolution process has 
caused the origin to be shifted by 3 grid units both in the 
X and Y direction. 

Figure 7 gives the gravity anomaly of an isolated 
prismatic model with its position indicated by dashed lines. 
The model has a length of -1.0, a width of 14. 0 and a 
thickness of 1.0. The upper surface is located at a depth 
of 3.0. Although the anomaly inflection line roughly 
corresponds to the vertical boundaries of the prism, the 
third vertical derivative of the anomaly enhances this 
information considerably (Figure 8). The maxima and minima 
flank the zero contour that marks the edges ot the body with 
reasonable accuracy. Unlike in magnetics, where there is a 
preferential shift of the contours towards one of the poles, 
depending on the direction of magnetization, the zero 
contour for gravity data shows symmetry around the body. 

The anomaly of two prisms in contact is shown in Figure 
9. Dashed lines indicate the position of the models. They 
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Figure 6. The third vertical derivative of a segment of a gravity 
anomaly generated for a three-dimensional rectangular prism (contour 
interval in mgals/di 3 ) 
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Figure 7. The gravity anomaly for a 
prism (contour intervals = 3 mgal) 


three-dimensional rectangular 
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Figure 8. The third vertical derivative of the gravity field for 
a three-dimensional rectangular prism (contour interval-120 mgal/di 1 ) 
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Figure 9. The superimposed gravity anomaly for two three-dimensional 
rectangular prisms in contact (contour interval ■ 10 mgal). 
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both have a depth of 3.0 and their density contrast is 0.2 
gn/cn . The map of their superimposed anomaly (Figure 9) 
gives the impression of being due to a single source with no 
apparent indication of the inner boundary. On the other 
hand, the map of the third vertical derivative of the field 
(Figure 10) reveals the existence of two sources. The zero 
contour closely follows the outer boundaries with reasonable 
accuracy, and the contact region is delineated by a well 
defined low trend. The absence of the zero contour along 
the contact is caused by the interference of the extrema 
from each model. This is a persistent problem also in 
magnetics for cases in which the distance between the bodies 
is approximately equal to or less than the depth of burial 
(Jain, 1974) . 

To simulate an area with several sources, a model with 
rectangular prisms of various orientations and sizes was 
evaluated. All the bodies have thicknesses of 1.0. Their 
upper surfaces are located at a depth of 1.5. Figure 11 
gives the superimposed anomaly of the different sources. 
The relative positions of the bodies are indicated by the 
dashed lines. The existence of the bodies is easily 
inferred from the map of the anomaly. However, the exact 
locations of their boundaries are not quite apparent from 
the map. On the other hand, in the map of the third 
vertical derivative (evaluated in the wavenumber domain), 


the sources are marked more accurately 


with the zero 
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Figure 10. The third vertical derivative of the superimposed 
anomaly for two rectangular prisms in contact ( contour 
interval * 10 mgal/di 3 ). 
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Figure 11 The superimposed gravity anomaly for five prismatic 
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Figure 12. The third * ical derivative of the 

superimposed anomaly for five prismatic sources with 

different sizes and orientations, (contour interval=10 mgal/di ) 
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contour( Figure 12). It is also noted that the orientation 
and relative positions of the various sources are well 
preserved by the interpre tat ional technique. 

Pete rmination of Depth to the Source 

It has been shown in Chapter II that the third vertical 
derivative may be employed to obtain the depth of a 
prismatic source. Equation 1 is, in fact, in a form which 
would allow one to estimate the depth to the top of an 
Isolated source. However, it has also been noted that, 
unlike in magnetics, the continuous distribution of density 
with depth requires the vertical extent be restricted to 
finite values. Thus, the implications of this on the value 
of the factor, C (Equation 1), has been studied in detail 
for two-dimensional prismatic models. 

The general procedure followed in this study involved 
measuring the horizontal distance between the maxima and 
minima of the third vertical derivative of a gravity field 
(see Figure 3) and then studying its relationship with the 
depth to the top of the source. This was systematically 
carried out by employing numerous models with various 
widths, thicknesses and depths. 

A sample of the depth estimations for two-dimensional 
models is given in Table 4. The thickness varies from 0.2 
to 12.0, while the depth and width are kept constant at 2.0 
and 12.0, respectively. 


For very thick sources (greater 



TABLE 4 


EFFECT OF THICKNESS ON DEPTH ESTIMATION 


Thickness 
( sampling unit ) 

0.2 

0.4 

0.6 

0.8 

1.0 

1.2 

1.6 

2.0 

4.0 

8.0 
12.0 


Distance between 
maxima and minima 
( sampling uni t ) 

1.8 

1.8 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.4 

2.5 
2.5 


than 8.0) the factor C in Equation 1 is 0.80. Fo*- bodies 
with intermediate thickness (0.6 to 2.0) it is 1.00. As the 
thickness decreases (below 0.6), C increases to 1.14. The 
general trend in the value of C is to be expected, since for 
decreasing thicknesses, the density distribution slowly 
approaches a distribution equivalent to magnetic monopoles 
at the top surface of a bottomless prismatic model (Chapter 
II). For practical purposes, a value of C equal to 1.07 
will, in general, give depth estimates with error bounds of 
less than IX provided the source is not thicker than two 
sampling intervals. Similar results are obtained when the 
above procedure was repeated for various values of the 
width. Some of the results are shown in Table 5. 
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TABLE 5 

EFFECT OF WIDTH ON DEPTH ESTIMATION 


Width 

(sampl ng unit) 


Distance between 
maxima and minima 
( sampling unit ) 


4.0 

6.0 
12.0 
18.0 
24.0 


2.0 

1.8 

1.8 

2.0 

2.0 


TABLE 6 

DEPTH ESTIMATIONS FOR A FIXED SAMPLING INTERVAL 


depth 

( sampling unit ) 

0.0 

1.0 

2.0 

3.0 

4.0 

5.0 

6.0 


Distance between 
maxima and minima 
( sampling unit ) 

0.4* 

1.2 

2.0 

3.0 

4.0 

5.1 
6.3 


*curve very oscillatory 


The range of depths the interpretatlonal technique will 
resolve for a pre-set sampling Interval was examined. Table 
6 gives selected results for two-dimensional models over a 
range of depths extending from 0.0 to 6.0. The thickness 
and width are kept fixed at 1.0 and 12.0, respectively. If 
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the value of C assumed earlier (1.07) is used, the depth 
estimates for models located between depths of 2.0 and 5.0 
are accurate to within about 9%. However, very shallow 
sources (less than 1.0) suffer from large errors primarily 
indicating that the selected grid interval was too large for 
such shallow depths. 

The results given in Tables 4, 5, and 6 cover a wide 
range of physical parameters that are of practical interest. 
If, for example, in a field survey the sampling interval is 
set to 10 km then the models cover depths extending from 0 
km to 60 km. The widths considered extend from 40 km to 240 
km. For most cases, the percentage error is less than 9%. 

It is important to realize that the above results were 
carried out for isolated sources. In cases where continuous 
sources (see Figure 12) are encountered, the interference of 
the maxima and minima of the various sources is to be 
expected. If it is impossible to select a representative 
isolated anomaly in such a region of interest, avoiding 
extrema on the side closer to the boundary of a major 
source would improve the reliability of the results. 
Moreover, as has been mentioned earlier, it is imperative to 
remove the undesirable effects of small, shallow sources 
prior to applying the operator. 


CHAPTER IV 


ANALYSIS OF 1° BY 1° MEAN FREE-AIR GRAVITY DATA 

In recent years, the theory of global tectonics has 
greatly improved our understanding of the dynamic processes 
at lithospheric plate boundaries which are believed to be 
formed by zones of accretion, subduction, and faulting. 
These zones, in turn, correspond to rift systems, trenches 
and transform and/or transcurrent fault zones, respectively 
(Wilcox, 1974). Linear zones of intense seismic activity 
and magnetic anomalies have, in the past, been used to 
identify and delineate such boundaries (Isaacs and Oliver, 
1968; Pitman, 1968). More recently, Wilcox (1974) and Kaula 
(1969, 1971) have utilized global free-air gravity maps for 
the same purpose. 

By comparison, little is yet known about the current 
tectonic regimes within the interior of continental plates. 
Earthquake seismology has not been fully utilized for 
detailed studies of parameters such as focal depth, state of 
stress, stress drop, and fault plane solutions, partly 
because only a small percentage of the world's total 
earthquakes contribute to the intraplate seismic activities 
(Sykes, 1978). Moreover, the seismograph station coverage, 
in most areas, is not adequate for detailed studies. This 
was especially true in the eastern and central United 
States, where earthquake stations were separated by a 
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typical distance of about 300 km until recently (Sykes, 

1978) . Hence, the examination of historical records and 
present-day seismic activities in this region has, for 
example, lead to the emphasis of certain trends of 
seismicity over others (Sbar and Sykes, 1973; Bollinger, 
1973, 1975, Woollard, 1958; Richter, 1959; King, 1970). 

Thus, at present, there is the need to supplement the 
investigations of seismic patterns with the study of 
lineaments revealed by satellite imagery, magnetic and 
gravity data, and various types of geologic information. In 
this chapter we will attempt to show how the third vertical 
derivative of the 1° by 1° mean free-air gravity data may be 
utilized to enhance certain linear features that may have a 
bearing on the identification of ancient intraplate 
boundaries, fault zones, suture zones, continental rift 
arms, or other major tectonic features. 

Global and regional mean free-air maps have over the 
years provided significant assistance in the examination of 
tectonic processes both at lithospheric boundaries and in 
their interiors. One of the major assumptions in such 
studies has been that sufficiently large areas are in 
isostatic equilibrium resulting in the mean free-air gravity 
value being approximately equal to zero (McGinnis, et al., 

1979) . Locally, however, it may be different from zero. 
This was inferred from the analysis of gravity data, and 
from various quantitative and qualitative analyses of 
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lithospheric bending under pressure (Gunn, 1949; Helskanen 
and Melnesz, 1958; McGinnis, 1970). It has been shown that 
the earth's crust behaves like an elastic plate under a 
topographic "load" and requires a relaxation distance of 
about 125 km to 200 km , depending on the geophysical 
parameters assumed (McGinnis, et at- 1979). 

Kaula (1969, 1972) described the tectonic classification 
of the global gravitational field on the basis of the 
magnitude and extent of mean free-air gravity anomalies for 
5° by 5° areas primarily derived from satellite data. He 
identified nineteen areas on the earth mostly over trenches, 
island arcs and active oceanic rises, as mainly positive, 
and fourteen areas primarily over oceanic and continental 
basins, as markedly negative. He also determined that the 
mean free-air gravity anomaly over the eastern United States 
is zero. Wilcox (1974) in a similar study utilized 1° by 1° 
mean free-air gravity values to identify patterns that are 
indicative of plate boundaries. He showed that a narrow 
belt of intensely negative values usually follows the 
present-day subduction zones (or trenches), while a belt of 
intensely positive values usually follows ocean rises and 
rift systems. 

Specifically, within the continental interior of the 
United States several attempts have been made to study the 
implications of regional iree-air gravity anomaly (Woollard, 
1976; McGinnis, 1970). A transcontinental gravity anomaly 
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profile extending from New York to Los Angeles showed that 
the United States is in lsostatic equilibrium (Woollard, 
1958). Similar results were also reported by Heiskanen and 
Melnesz (1958). In fact, further studies by Woollard (1962) 
and McGinnis (1970) have shown that the relaxation distance 
for the continental interior is about 250 km, and the 
average free-air anomaly value over such large areas would 
be expected to be approximately equal to zero. Hence, the 
deviation of the mean values from zero indicates isostatic 
imbalance. The resulting lineations on the anomaly map may 
then be utilized to delineate intraplate boundaries or to 
outline major crustal blocks. 

In the next sections, some major lineaments seen in the 
most recent 1° by 1° mean free-air gravity data of eastern 
and central United States will be presented. The 
interpretational procedure will be enhanced by applying the 
vertical derivative filter developed in the previous 
chapters. Such linear features, apart from their tectonic 
significance are also increasingly being recognized as 
possible indicators of mineral deposits and sites of 
hydrocarbon accumulations (Saunders, et al., 1976). 
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Correlations of the Mean Fr ee-Al r Gravity Da ta 

with Structure 

The gravity field data used in this paper are part of the 
most recent 1° by 1° mean free-air gravity anomalies of the 
earth compiled by Rapp (1978). The terresterial gravity 
values for the eastern and central United States are the 
updated and expanded versions of the 1978 data set of the 
Defense Mapping Agency Aerospace Center (DMAAC) and the 
Lamont-Doher ty Geological Observatory. Over some of the 
continental margins and the ocean areas the data are 
obtained from GEOS-3 satellite altimetry. The average 
global accuracy of Rapp's (1978) data set is +15mgal; 
however, he noted that in the eastern and southeastern 
United States, where gravity data are dense, the accuracy is 
expected to be much better (Rapp, 1978). 

To facilitate computer plotting and contouring, the data 
were interpolated onto a square grid of about 1/4° by 1/4° 
using a program based on the cubic-spline algorithm 
published by Davis (1972). The output was then 
automatically contoured at a 10 mgal interval using the 
Alber's equal area projection to produce the free-air 
anomaly map (Figure 13). 

The overall pattern of the 1° by 1° free-air anomalies 
seems to follow the major tectonic provinces, especially in 
the Appalachian foldbelt. The map is very similar to the 1 0 
by 1° mean free-air map of Woollard (1976) and to the map 



Scale 1:15,000,000 


figure 13. 1 by 1 free-air gravity data for the central 

and eastern United States. See the text for the meaning 
of the symbols (contour interval=10 ragal). 
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constructed by McGinnis, et al., (1979) for values averaged 
over twelve quadrangles (1° by 0.8°). The anomaly values 
range from over +30 mgal in the foldbelt to below -50 mgal 
over the continental margins. 

For ease of discussion, the prominent highs and lows that 
occur over major tectonic regions (Figure 14) are marked 
with letters on Figure 13. In the eastern central region 
linear belts of negative and positive anomalies trend in 
the northeastern direction. The inner Piedmont (AD) and the 
deepest parts of the Appalachians are associated with a 
mildly negative anomaly. On the other hand, a strongly 
positive gravity high follows the east central Piedmont 
across North Carolina and Virginia (CP). Another prominent 
high over the Blue Ridge Mountains (BR) is found to the west 
on line with the Kane gravity high (KG), in Pennsylvania, to 
the north. The southeastern extension of the large gravity 
high passes through Georgia (GF) and loops around Florida 
and heads westward along the Gulf Coast. To the north, a 
small remnant of the Scranton gravity high joins the 
positive gravity anomaly associated with the Adirondack 
Uplift (SA) resulting in a change of the general trend to 
the north. Over the coastal areas and continental margins 
to the east, the anomalies are predominantly negative (CM). 

The belt of large, positive gravity anomalies in the 
northwestern corner of Figure 13, is a segment of the so- 
called aidcontinent gravity high (MC), which is even more 
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Figure 14. Location of major tectonic features in the 
central and eastern United States. 1. Piedmont; 2. Brevard 
Fault Zone; 3. Blue Ridge;4. Adirondack Uplift;5. Gettysburg 
Basin;6. Newark Basin;7. Connecticut Basin;8. Michigan Basin; 
9. Wisconsin Arch;10.Cincinati Arch; 11. Illinois Basin ;12. 
Mississippi Embavment;13. Atlantic Coastal Plain;14. Blake 
Fracture Zone;17. Ouachita Mountains. Modified after 
McGinnis, et al . (1979). 
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prominently seen in the Bouguer anomaly map of the United 
States (Woollard, 1964). It is flanked by a negative 
anomaly over the Wisconsin Arch (WA) and by a positive 
anomaly over the Michigan Basin (MB). The Mississippi 
Embayment (ME), to the south, is associated with a broad, 
mildly positive anomaly that joins the Kentucky gravity high 
(KH) to the east. There is a gravity low over the Ouachita 
foldbelts to the west (QF). A negative anomaly corresponds 
to the Cincinnati Arch (CN). 

In general, the association of regional changes of 1° by 
1° free-air anomaly with major tectonic regimes seems 
apparent. As has been noted by King (1970) and Woollard 
(1976), these anomalies should largely reflect the property 
of the underlying masses, since the surface relief in the 
central and eastern United States is, on the average, 
moderate. There are some notable exceptions to this over 
areas with high elevation, such as the small gravity high to 
the southeast of the Kentucky gravity high. One also notes 
that positive anomalies occur over diverse features such as 
the Michigan Basin and the Adirondack Uplift, while negative 
anomalies are observed over the Wisconsin Arch and the 
Cincinnati Arch. This probably Indicates that the anomalies 
are not caused by shallow sources. 

There have been some attempts in the past to understand 
the nature c * the sources that contribute to regional 
gravity anomalies derived from ground measurements or from 
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Che use of satellite orbital data. Woollard and Khan (1970) 
conducted a conventional half-width analysis on 12° by 12° 
free-air gravity data from satellite measurements, and 
showed that for the major anomalies considered, the centers 
of the contributing masses have maximum depths that are less 
than 1000 kms in all cases, and not more than 150 kms in 
most cases. Woollard (1976) has also made a comparative 

O 0 

study of 1 by 1 free-air gravity anomaly with the change 
of geophysical parameters along the 45°N parallel across the 
United States. He showed that the pattern of changes of the 
anomaly arj related to the integrated effect of the mean 
crustal velocity, the values of the crustal thickness, and 
the upper mantle velocity. In regions where these 
parameters are high, there are positive gravity anomalies. 
On the other hand, where they are low, negative anomalies 
are observed. Hence, it may be concluded that the 1 ° by 1° 
free-air anomalies mainly manifest the overall long 
wavelength changes in crust and upper mantle structure. 



53 


Correlation of Data with Regional Lineaments 

It has been shown earlier that the zero contour of the 
third vertical derivative closely outlines the boundaries of 
an isolated prismatic source. For a continuous distribution 
of models, the areas of contact were shown to be delineated 
by well-defined low or high trends. It has also been noted 
that some anomalies in the mean free-air gravity map (Figure 
13) are associated with major tectonic regimes (Figure 14). 
In fact, one also observes that some linear features in the 
anomaly map are defined by the zero contour together with 
some major trends of highs or lows. In this section, the 
vertical derivative operator will be employed to enhance 
such linear features in the 1° by 1° mean free-air gravity 
data for the central and eastern United States. 

Sampling theory implies that, the scale of features for 1° 
by 1° (or about 110 km by 110 km) data is roughly 220 km. 
However, the gravity anomaly map (Figure 13) and its 
amplitude spectrum (Figure 15) indicate the presence of some 
higher wavenumber features which might lead to instability 
problems when obtaining the third vertical derivative. Such 
problems find their cause in the tendency of the vertical 
derivative operator to amplify the high wavenumber content 
of che field spectrum, which is caused by small, shallow 
sources or by random errors introduced during the data 
acquisition, digitization, and/or interpolation processes. 
A detailed Inspection of the amplitude spectrum showed that 
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Figure 15. The amplitude spectrum of the gravity data for the central and 
eastern United States (contour interval = 5 units). 
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removing wavelengths shorter than about 250 km (or .04 
cyclea/dl) was adequate for our signal processing schemes. 
The data were convolved with a set of 13 by 13 low-pass 
filter coefficients that were obtained using an algorithm 
given by Lavin and Devane (1970). The normalized amplitude 
response for the filter (Figure 16) shows that it is 
circularly symmetric for wavenumbers less than about 0.1 
cycles/dl. Hence, it should not introduce any directional 
effects in the output. 

The filtered version of the free-air data was convolved 
with the 7 by 7 set of coefficients of the third vertical 
derivative operator that was developed in the previous 
chapter (Figure 17). We recall that the operator gives the 
maxima and minima of the anomaly-curvatures near the 
boundaries of the isolated body. The distance between the 
extrema was related to the depth of the top surface of the 
source, while the zero-contour was found to delineate the 
outlines of the body. However, in view of the fact that the 
observed gravity field of the eastern and central United 
States is the result of the superposition of highly complex 
mass distributions, the position of the maxima and minima 
from each source cannot in most instances be accurately 
isolated. In fact, the interference between the extrema 
will have the tendency to complicate the already complex 
picture, as is well illustrated in Figure 17. Hence, the 
attempt to interpret the results will be directed only 
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Figure 16. The normalized amplitude response for the low-pass 
filter (contour interval = 3db). 
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Cowards the utilization of the zero-contour of the third 
vertical derivative of the data to identify linear features 
in the area. These lineaments may then outline the position 
and areal extent of major crustal blocks. In order to 
facilitate this interpretation the extrema of the vertical 
derivative of the data were removed through low-pass 
filtering, and the resulting map is shown in Figure 18. 

It is noted that the maps shown in Figures 13, 17, and 18 
have become progressively smaller, though the scale is kept 
constant. This behavior is caused by the loss of data 
around the edges during the convolution process. 

It should be emphasized that the filtering procedure used 
does not necessarily eliminate the effect of shallow 
features, since any particular range cf wavenumbers of the 
spectrum contains the contributions from various sources at 
different depths. This implies that the gravity effect of 
regional structures in the sedimentary layer is superimposed 
on the effects of deeper structures. Likewise, the nature 
of the 1° by 1° free-air gravity data poses some problems. 
Even though the data set is the most accurate to date, there 
are still some problems in the data acquisition and 
averaging processes that are being worked out (Rapp, 1978). 

For this reasons, the analysis, in this section, will be 
solely aimed at demonstrating the advantages and limitations 
of the third vertical derivative operator for such a data 
set. The overall results should, thus, be taken as 
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Figure 18. The smoothed version of the third vertical 
derivative of the gravity data for the central and eastern 
United States (contour interval = 50 mgal/di 3 ) See the 
text for the meaning of the symbols. 
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Prominent linear features that are enhanced by the 
vertical derivative filter are shown by solid lines in 
Figure 18. These lineaments are inferred from zones 
delineated by the zero contour, and from areas along which 
high or low trends are aligned or terminated. Only the very 
distinct features that appear to have possible bearing on 
identifying major crustal blocks are marked. The overall 
picture in Figure 18 seems to show a bias toward the 
dominance of major lineaments in the northeast-southwest 
directions. Affleck (1963, 1970) reported that the 
aeromagnetic data showed a similar trend as well as a few 
others. (He also reported that the vertical derivative 
filter was employed to enhance the interpretation of the 
data ) . 

In order to form some framework for the analysis we have 
assigned a letter to each of the linear features (Figure 
18). Lineament A, which is, perhaps, the most persistent 
one, correlates with the so-called New York-Alabama magnetic 
lineament (King and Zeitz, 1978) all the way from the 
southern part of New York to eastern Tennessee. It swings 
to the southeast in Alabama, where the "magnetic lineament 
was absent and was actually augmented using the Bouguer 
gravity anomaly (King and Zeitz, 1978). We also note that 
lineament A tends to follow the west side of the Appalachian 
gravity low (AD in Figure 13) and also separates a zone of 
predominately north trending gravity anomalies to the 
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northwest from a zone of northeast trending anomalies to the 
southeast. This Is clearly seen in Figure 17. Moreover, as 
has been noted by Watkins (1964) and King and Zeitz (1978) 
for the New York-Alabama lineament, A seems to bound an 
active seismic zone over the eastern part of the Appalachian 
basin extending from New England to eastern Tennessee 
(Figure 19). It also parallels the New Madrid earthquake 
zone to the T est. King and Zeitz (1978) have suggested that 
the zone may represent an ancient intraplate boundary 
between the highly deformed Appalachian foldbelt and the 
stable craton of the continental interior. 

Lineament B traces the Brevard fault zone throughout its 
mapped lergth that extends from southern Virginia to the 
coastal plains of Alabama .(Reed and Bryant, 1961; King, 
1977). Saunders and Hicks (1976) have noted that satellite 
imagery suggests furtner extensions northward into Maine and 
southward to the Gulf coast. However, lineament B is 
truncated in Virginia and merges with New York-Alabama 
lineament to the south. 

Along the eastern boundary of the Piedmont ("the fall 
line"), there is a belt of earthquakes that extends roughly 
from New Jersey to Alabama, along which many major faults 
have been recognized (Agarwal and Sykes, 1978; Higgins et 
al., 1974). A focal mechanism solution for an earthquake 
around Wilmington, Delaware, showed dip-slip movement on a 
nearly vertical fault oriented in a northeast direction 
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(Sbar, et al., 1975). Higgins ,et al. (1974) have also 
noted that the pattern of aeromagnetic anomalies suggest a 
system of faults that runs from New York to Virginia. 
Lineament C appears to correlate with such a zone of 
faulting and seismic activity. In the north this linear 
feature intersects lineament D, which traces the axis 
defined by the Appalachian Salient between South Mountain 
and the Newport gravity and magnetic high. Further 
extension of lineament D to the northwest follows the 
Tyrone-Mt. Union lineament. 

In the northwestern part of the region, lineament G 
appears to bound the New Madrid earthquake zone. This 
seismic region is the most active in the central and eastern 
United States. Historic records that include two large 
earthquakes of 1811 and 1812, and recent activities are 
concentrated along the northern part of the Mississippi 
Embayment, which is associated with northeast trending 
faults, and large magnetic anomalies ( Syke s , 1 9 7 8 ) . 

Lineament "H," in the south, partly defines the Blake 
Fracture Zone, an ocean floor topographic feature where a 
left lateral offset (100km) of magnetic lineaments has been 
observed (Pitman and Taiwan i , 1 9 7 2 ) . On land, extensions of 
this linear feature are inferred from satellite Imagery and 
geophysical data. Bollinger ( 1972 , 1973 ), ir. a study of 
historic earthquakes and recent activities defines a South 
Carolina-Georgia zone of northwesterly seismic trend which 
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is not far from an extension of the fracture zone. The 
trend is also delineated by several Igneous intrusions 
(King, 1961; De Boer, 1967), and by some prominent magnetic 
anomalies indicative of mafic and ultramafic rocks at depth 
(Popenoe and Zeitz, 1977). It is also interesting to note 
that the seismic area close to Charleston, South Carolina, 
falls close to the intersection of this linear feature with 
lineament C. 

At this stage of the project, one may attempt to identify 
the crustal blocks that are bounded by the major lineaments 
keeping in mind the problems brought on by the superposition 
of shallow and deep sources. For example, King and Zeitz 
(1978) have recognized the aseismic area northeast of the 
Mississippi Embayment as a more stable crustal block. It is 
bounded by lineament A in the east and the New Madrid 
s°ismic zone in the west. The southern boundary is the zero 
contour of the Mississippi Embayment. Its border in the 
north is not clear from the map though an extension of 
lineament D is a possible candidate. Another crustal block 
is, perhaps, what Sykes (1978) calls the Blue Ridge Dome. 
This block is roughly bounded by the zero contour (Figure 
18) marked by lineaments A, B, C, D, and an approximate 
ex nsion of lineament H on the south. It is possible that 
one will be tempted to continue such a scheme over the whole 
region. For instance, the shaded areas in Figure 18 that 
are delineated by the zero contours are some likely 



candidates for crustal blocks. One may then attempt to 
study the seismic, magnetic and other data for the 
indicated blocks to determine the difference between their 
average properties. However, as it has been stressed 
earlier, any detailed and reliable analysis should wait 
until the influence of the sedimentary layer is identified 


and eliminated 



CHAPTER V 


CONCLUSIONS 

The major intent of this study has been to investigate 
the conditions whereby an inte rpre ta t ional technique that 
has been extensively used for aeromagnetic data analysis may 
be extended to gravity field Interpretation. The 
technique was based on the fact that, for an isolated, 
bottemless rectangular prism, the zero contour of the second 
vertical derivative of its magnetic field closely outlines 
its top surface (Vacquier et al., 1958). Moreover, the 

horizontal distance between the plotted maxima and minima of 
the derivative gives an estimate of the depth to the source 
(Jain, et al., 1974). An examination of Poisson's theorem 
indicated that an extension of the method to gravity field 
analysis was possible. It implied that the second vertical 
derivative of the magnetic field is equivalent to the third 
vertical derivative of the gravity field. However, The 
continuous distribution of density with depth required that, 
when using the gravity field, the prismatic model had to be 
restricted to a finite thickness in order to avoid a 
singularity. Hence, a detailed model study was carried out 
to determine the effects of the thickness factor on the 
in terp r e ta t iona 1 technique. 

The rectangular prismatic model with finite thickness was 
examined extensively because of its wide applications in the 
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Investigation of varied geological structures. The space 
domain expression of its gravity effect showed that its 
parameters were contained in coupled terms . This made the 
analytical evaluations of the third vertical derivative 
rather inter ac table . Hence, for the sake of simplicity, 
most of the analysis was carried out in the wavenumber 
domain. In some cases, it was necessary to employ a 
derivative operator designed in the spatial domain in order 
to circumvent severe truncation effects. 

An examination of the sensitivity of the interpre tat ional 
technique to sampling interval demonstrated that the finer 
the grid spacing the more accurate the computed results. 
However, a sampling internal of about 0.2 times the width 
for the numerical models gave adequate results in outlining 
the top surface of the source, while a value of 0.1 was 
required for estimating their depth of burial. In cases 
where the grid spacing was coarser, interpolation to finer 
spacing using the bi-cubic spline interpolation in the 
spatial domain and padding with zeroes in the wavenumber 
domain were found effective in improving the results. 

The model study covered a wide range of physical 
parameters that were of practical Interest. If, for 
Instance, a grid interval in a field survey is assumed to be 
10 km, the range of thicknesses studied varied from 2 km to 
200 km, the depths considered extended from zero to 70 km, 
and the width varied from 20 km to 240 km. The overall 
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results demonstrated that, for thicknesses ranging from 0.2 
to 20.0 sampling Intervals, error bounds of less than 8 
percent are to be expected on outlining the edges of the 
source. On the other hand, error bounds for the depth 
estimate would be less than 9 percent If the thickness Is 
smaller than twice the grid interval. 

The versatility of the lnterpre tational techniqu* was 
illustrated by evaluating the superimposed anomalies oi 
several three-dimensional rectangular prisms. It was 
demonstrated that the zero contour of the third vertical 
derivative of the gravity field closely outlines the top 
surface for isolated sources. In cases where the models 
were in contact, the boundaries were delineated by well 
defined low or high trends. The orientations and relative 
locations of the various models were also well preserved by 
the technique. However, the depth estimation for sources in 
contact was made difficult by the interference of their 
extrema. In such cases, using extrema that are away from 
the contact region gave bet er results. 

Finally, an attempt was made to analyze the most recent 
1° by 1° free-air gravity data for the central and eastern 
United States. The anomaly map produced was found to be 
very similar to that of Woolla rd ( 1 97 6 } and McGinnis, et al., 
(1979). The general pattern of the anomalies appears to 
follow the major tectonic regimes. It was understood that 
the major high and lows primarily reflect the effect of the 
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underlying masses, since Che mean surface elevation in Che 
region is moderaCe. Moreover, Che observation of, for 
instance, large positive anomalies over both basins aril 
uplifts (e.g; che Michigan Basin and Adirondack. Uplift) 
seems Co indicate chac Che sources of me anomalies are not 
shallow in nature. In fact, Che overall picture .appears Co 
support Che suggestion by Woollard (1976) thaC Che free-air 
gravity anomaly reflects the integrated effect of the mean 
crustal and upper mantle structure. 

The inte rpre tat ional technique developed in this pa, sr 
was employed to enhance the linear features observed in the 
1° by 1° free-air gravity data. Prior to this, the high 
wavenumber content of the field spectrum was removed through 
low-pass filtering. However, it is likely that the 
filtering procedure did not necessarily eliminate the effect 
of all shallow sources. Hence, only a preliminary 
analysis of the most distinct lineaments was attempted. In 
spite of this res t r let ion , the correlation of the inferred 
linear features with some seismic trends and magnetic 
lineaments was demonstrated. Some possible crustal blocks 
were also delineated by the zero contour of the third 
vertical derivative of the field data. 

Further werk should be done on the relationship between 
the Inferred crustal blocks and the basement structure in 
the region. This probably involves isolating and removing 
the effect of the regional stuctures in the sedimentary 
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layer. Further work should also be done in comparing the 
seismic, magnetic and other geophysical data for the 


inferred 


crustal blocks. 
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